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Abstract: We present a detailed study of the charged current Drell-Yan process, which 
includes the exact 0(a) electroweak corrections properly matched with leading- log effects 
due to multiple-photon emission, as required by the experiments at the Tevatron and 
the LHC. Numerical results for the relevant observables of single W boson production at 
hadron colliders are presented. The impact of the radiative corrections and of some sources 
of theoretical uncertainty is discussed in detail. The calculation has been implemented in 
the new version of the event generator HORACE, which is available for precision simulations 
of the charged current Drell-Yan process. 
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1. Introduction 

At hadron colliders, such as the Fermilab Tevatron and the CERN LHC, the production 
of a high transverse momentum lepton pair, known as Drell-Yan process plays an 
important role: it allows, in the charged current channel, a high precision determination of 
two fundamental parameters of the Standard Model, namely the mass and the decay width 
of the W boson Q; it provides, both in neutral and charged current channels, stringent 
constraints on the density functions which describe the partonic content of the proton Q ; 
it can be used as a standard reference process and therefore as a luminosity monitor of the 
collider [[|, M. Furthermore, it represents a background to the search for new heavy gauge 
bosons ||. 

The accuracy in the determination of the theoretical cross section has greatly increased 
over the years. The calculation of next-to-leading order (NLO) QCD corrections |7| has 
been one of the first test grounds of perturbative QCD. Next-to-next-to-leading order 
(NNLO) QCD corrections to the total cross section have been computed in ref. [||, but 
differential distributions with the same accuracy have been obtained only recently in ref. || . 
The size of the NNLO QCD corrections and the improved stability of the results against 
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changes of the renormalization/factorization scales raises the question of the relevance of 
the 0(a) electroweak (EW) radiative corrections, which were computed, in the charged 
current channel, first in the pole approximation [10, 11] and then fully in refs. [12, 13, 14, 



1§. 

A realistic phenomenological study and the data analysis require the inclusion of the 
relevant radiative corrections and their implementation into Monte Carlo event generators, 
in order to simulate all the experimental cuts and to allow, for instance, an accurate 
determination of the detector acceptances. The Drell-Yan processes are included in the 
standard QCD Parton Shower generators HERWIG and PYTHIA [16, 17]. Recently there 
have been important progresses to improve the QCD radiation description to NLO, which 
has been implemented in the code MC@NLD [l~£[| . Another important issue is the good 
description of the intrinsic transverse momentum of the gauge boson, which can be obtained 
by resumming up to all orders the contributions of the form a s \og(p±/m w ). The generator 
RESBDS fl~9[| , used for data analysis at Tevatron, includes these effects. 

If the inclusion of QCD radiation is mandatory for the simulation of any process at 
a hadron collider, one should not neglect the impact of EW corrections on the precision 
measurement of some Standard Model (SM) observables, like the W boson mass and decay 
width. For instance, the generator WGRAD [14] includes the exact 0(a) EW corrections, 
which have been shown to induce a shift on the value of m w extracted from the Tevatron 
data of about 160 MeV in the muon channel [|| , mostly due to final-state QED radiation. 
In view of the very high experimental precision foreseen at the LHC (Am w ~ 15 MeV), 
final-state higher-order (beyond 0(a)) QED corrections may induce a significant shift, as 
shown in ref. p^] . Some event generators can account also for multiple-photon radiation: in 
the published version of HORACE [H, H final-state QED radiation was simulated by means 
of a QED Parton Shower [22]; the generator WINHAC [p3j | uses the Yennie-Frautchi-Suura 
f24|1 formalism to exponentiate final-state- like EW 0(a) corrections; finally, the standard 
tool PHOTOS p|] can be used to describe QED radiation in the W decay. 

A first attempt to study the combined effect of EW and QCD corrections has been 
presented in ref. pGfl . 

The predictions of WINHAC and HORACE and of WINHAC and PHOTOS, for TV-decay, have 
been compared in the papers of ref. [27|. A detailed series of tuned comparisons between 
different EW Monte Carlo generators have been done in ref. [28], in order to check the 
reliability of different numerical predictions, with 0{a) accuracy, in a given setup of input 
parameters and cuts. 

Since the Drell-Yan events can be used, in principle, to determine the collider lumi- 
nosity at a few per cent level, the theoretical cross section must be known with the same 
accuracy, requiring also the inclusion of 0{a) EW corrections. Furthermore, the 0{a) EW 
contributions give large corrections to the tails of the transverse mass and lepton tranverse 
momentum distributions, because of the presence of large EW Sudakov logarithms [13, 14]. 
These regions are important for the search of new heavy gauge bosons. 

The aim of this paper is to present a precision EW calculation of the charged current 
Drell-Yan process, which includes the exact 0(a) EW matrix elements properly matched 
with leading-logarithmic higher-order QED corrections in the Parton Shower approach. 
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The matching here presented of perturbative corrections with Parton Shower, which is a 
topic of great interest in modern QCD simulations (^9|, is the first example of such an 
application in the EW sector and is realized along the lines already presented in ref. pO|j , 
Several distributions of physical interest are analyzed, disentangling the effect of different 
classes of radiative corrections and discussing various sources of theoretical uncertainty. 
The calculation is implemented in the Monte Carlo event generator HORACE, which com- 
bines, in a unique tool, the good features of the QED Parton Shower approach with the 
additional effects present in the exact 0(a) EW calculation. This task is non trivial from 
several technical points of view and faces all the conceptual problems of developing a NLO 
event generator. 

The paper is organised as follows. In Section ^ we present the calculation of the 
0{ot) EW corrections to the partonic process ud — ► l + i>i(^). In Section [3| we describe 
the matching of the fixed order results with the QED Parton Shower. In Section ^ we 
present the computation of the hadron- level cross section o (^pp — > l + ui(nj) S j and discuss 
the subtraction of the initial-state collinear singularities to all orders. In Section |5] we 
present phenomenological results for several physical distributions and discuss the impact 
of EW O(a) and of higher-order QED corrections. Finally, in Section [6| we draw our 
conclusion and discuss possible developments of this work. 



2. Partonic process: matrix elements calculation 
2.1 Born approximation 

We consider the charged current Drell-Yan partonic process u{p\) d{p2) — > vi(ps) 1 + {pa)- 
This process is a weak charged current process and its amplitude, in unitary gauge, is 
proportional, at tree level, to the square of the SU(2)l coupling constant g: 
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where V u d is the CKM matrix-element, m w is the TV-boson mass and Y w is the W decay 
width, necessary to describe the W resonance region, s = {p\ +P2) 2 is the squared center- 
of-mass energy and = p^ +p%- The differential cross section, mediated over initial-state 
and summed over final-state spins, mediated over initial-state colours and in the limit of 
vanishing fermion masses, is given by 
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where u = {jp\ — Pa) 2 - The weak coupling g can be expressed in terms of a, the fine 
structure constant, and of sva.6 w , the sinus of the weak mixing angle, via the relation 
Ana = g 2 sin 2 9 W . The weak mixing angle is defined as cos6 w = m w /m z , where m z is the 
Z boson mass. 
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Figure 1: Some examples of one-loop virtual diagrams. 



2.2 The O(o) calculation 

The complete EW 0(a) corrections to the charged current Drell-Yan process have already 
been computed in refs. |l0|, 11, 12, [l^, 14, 15|. We have repeated independently the 



calculation and we summarize here its main features. 

The 0(a) corrections include the contribution of real and virtual corrections. The 
virtual corrections follow from the perturbative expansion of the 2 — ► 2 scattering amplitude 
M = M + M v a irt + ■■■ and contribute, at 0(a), with 2Re(M v J rt M* ) . The 0(a) virtual 
amplitude includes two contributions, namely the one- loop renormalization of the tree- 
level amplitude and the virtual one-loop diagrams. The real corrections are due to the 
emission of one extra real photon and represent the lowest order of the radiative process 
u(pi)d{jp2) — ► i / i(P3)l + (p4)'y(k). They can be further divided in soft and hard corrections, 
M\ = M s °^ + Aii ard . The former respect, by definition, the Born-like 2 — > 2 kinematics 
and can be factorized as |.M^°^*| 2 = ^sbI-MoI 2 where 5sb is a universal factor that depends 
only on the external particles. The total cross section includes soft and hard corrections 
and is independent of the cut-off used to define the two energy regions. Virtual and real soft 
corrections are separately divergent due to the emission of soft photons, but the divergence 
cancels in the sum of the two contributions. 



2.2.1 Virtual corrections 

The 0(a) virtual corrections to a 2 — * 2 reaction include contributions of counterterm, self- 
energy, vertex and box corrections. Few diagrams representative of the different kinds of 
corrections are depicted in figure and have been calculated using the packages FeynArts 
and FormCalc |3l|, [32j|. The numerical evaluation of the 1-loop integrals has been done 



using the package LoopTools2 [32], based on the library ff |33|]. We will write the 1- 
loop virtual amplitude as M v J rt = M% s + Ma lf + M v a ertex + M b ° x . We have explicitly 
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checked that our results numerically agree with those of refs. 15|; the fermion mass 



corrections, which are present in our calculation and neglected in refs. 15 1, turn out 
to be negligible from a numerical point of view. The mass of the fermions in the scalar 
1-loop integrals regularizes in a natural way the mass singularities due to the emission of 
a (virtual) collinear photon. The infrared divergence of the integrals has been regularized 
by means of a small photon mass A. 

The introduction of the W decay width in the propagator of the W boson is mandatory 
to describe the resonance region and to regularize the divergence due to the pole of the 
propagator. In order to account for the W width, we observe that the W propagator, in 
the 't Hooft-Feynman gauge, can be modified at 1-loop in the following way: 

HO- 



{-igT) (IWW + Sm 2 w + (s- m 2 w )5Z w ) (2.3) 

s - mfy + iT w m w ' s - 

where Hww(s) is the transverse part of the W self-energy corrections, Sm^ and 5Z W 
respectively the W mass and wave function renormalization constants. The two countert- 
erms cancel the divergences present in the self-energy corrections. We remark that the 
second factor l/(s — m 2 ^) in eq. (2.3) is not corrected by the decay width, to avoid double 
counting; we can check, by expanding the self-energy corrections about s = rri^, that the 
0(a) expression is regular for s — > m 2 ^. The contemporary presence of the decay width in 



the propagator and of the explicit 0(a) corrections does not yield any double counting |i3| : 
in fact the imaginary part of the self-energy corrections contributes only at the 2-loop level 
to the cross section and does not enter in the 1-loop virtual contribution 2Re (A4^ r *7Wo)- 
The W self-energy includes also contributions from diagrams having a W and a photon 
running in the loop, which develop a logarithmic divergence when s = m^- in the 1-loop 
virtual amplitude, the coefficient of log(s — m^) is gauge invariant; we can therefore reg- 
ularize the logarithmic divergence by replacing log(s — m^) — > log(s — m% + iT w m w ) in 
all the self-energy, vertex and box scalar integrals, without spoiling the gauge-invariance 
of the calculation. 

Examples of vertex and box corrections are depicted in figure |l[ The vertex diagrams 
with a trilinear gauge boson vertex and the box diagrams with a photon in the loop yield 
the logarithmic divergence at s = m^, previously discussed. The abelian vertex diagrams 
and the box diagrams with a photonic correction are infra-red divergent. All the vertex 
and box diagrams with a Z boson exchange yield the so-called EW Sudakov logarithms, 
namely terms like a log 2 (s/m|), whose importance grows for large invariant mass of the 
final state lepton pair, while they are almost negligible at the W resonance. 

The calculation has been repeated with two different gauge choices, namely the i?g 
with £ = 1 gauge and the background field gauge, with parameter q = 1 |34]]. The two 
results perfectly agree, and this is an important check on the calculation of the bosonic 
self-energy and of the non-abelian vertex corrections. 

Concerning the renormalization of the 1-loop amplitude, the UV divergences which 
appear from the virtual diagrams can be cancelled with the mass, Sm^, and wave func- 
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tion, 5Z W , renormalization constants of the W boson and by the renormalization of the 
two vertices Wud and Wl + vi. The latter include the charge renormalization and the wave 
function renormalization of the external fermions and of the W boson. Being the tree-level 
amplitude of eq. (2.1) proportional to g 2 , using the electric charge and the gauge boson 



masses (e, m w ,m z ) as input parameters, we write the bare coupling go = eo/s W) o and then 
replace it in terms of renormalized quantities and of counterterms go = e/s w (1 — <5e/e)/(l — 
Ss w /s w ). The electric charge counterterm is fixed by the request that in Thomson scatter- 
ing the renormalized charge is given by the fine structure constant; its expression depends 
on the quark masses running in the photon vacuum polarization, the value of which can 
be adjusted in order to make the running electric charge reproduce the value a(m? z ) p5| . 
The weak mixing angle counterterm is given as a combination of the mass counterterms of 
the W and Z bosons, following from its definition. The W boson mass and wave function 
renormalization constants are defined in the on-shell scheme. 

The choice of the input parameters of the SM lagrangian has an impact on the predic- 
tion of the physical observables. If ideally one were able to resum exactly the perturbative 
expansion, the predictions would be the same in any scheme. On the other hand the trun- 
cation of the perturbative expansion induces a dependence on the scheme, which is formally 
of higher order, but which can be numerically relevant. 

Some of the possible options in the gauge sector of the EW SM are: i) c^G^m^ 
ii) a,m w ,m z , hi) G^, m w , m z . The first option, used for LEP1 analyses, is based on the 
best measured EW quantities and minimizes the parametric dependence of the predictions 
on the inputs. In this input scheme, the value of m w is a predicted quantity. For W physics 
at hadron colliders, the other schemes are preferable because ra w is an input parameter. 

The second choice (a(O)-scheme) has the proper coupling a for the real photon emission 
diagrams and parametrizes the charged current coupling g as V 'Aira / 's w . 

The choice which seems to be more natural for W physics is the third one (G^ scheme), 
where the weak coupling g is related to the Fermi constant and to the W boson mass by 
the relation 

I = i4 (1 + Ar) (2 - 4) 

The quantity Ar represents all the radiative corrections to the muon-decay amplitude |36| . 
Introducing 

,2 a, fnr< ™2 



9 



4^26^774(1 - Ar) (2.5) 



in the tree level amplitude (2.1) we generate an additional contribution to the 0(a) correction 
proportional to Ar. Being the vertex between charged particles and photons proportional 
to gsin6 w , we can therefore introduce an effective electromagnetic coupling constant 



/ „ , V^Gfj, sin 2 Owin^ 
a G u = ( 2 -6) 



which is derived from eq. (2^) and is evaluated in tree- level approximation by setting 
Ar = 0. The effective coupling olq^ differs from the fine structure constant a, evaluated 
at zero momentum transfer, by higher-order effects. 
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(a) (b) (c) (d) 

Figure 2: 0(a) bremsstrahlung diagrams. 



2.2.2 Bremsstrahlung corrections 

The real radiative corrections to the charged Drell-Yan process, described by the amplitude 
Mi, are given by all the Feynman diagrams (figure ||) with the emission of one extra photon, 
from all the electrically charged legs of the Born diagram, including the internal W boson. 

The probability amplitude has been calculated in the unitary gauge with massive 
fermions. We integrate the squared matrix element over the whole photon phase space 
and split the allowed photon energy range in two intervals, [A, AS] and [AE, E max ]. The 
cut-off AE <C -^/i is chosen in such a way that the photon with smaller energy is considered 
soft and does not modify the 2 — > 2 kinematics of the Born amplitude. The small photon 
mass A has been introduced to regularize the infrared divergence. In this energy region the 
phase space integral, including the full angular integration, can be solved analytically. The 
result can be expressed in a factorized form, as 



f=u,d,e+ 



where the soft Bremsstrahlung factor, see e.g. |$7|], depends on the mass and electric charge 
of the external radiating particles and the phase-space region $7 is defined by the request 
that the photon energy Ey satisfies A < EL < AE. We have explicitly checked that the 
sum of the virtual and soft-real contributions is independent of the choice of the photon 
mass A, in the limit of small A values. 

In the hard energy region the phase space integration has been performed numerically, 
with Monte Carlo techniques improved by importance sampling to take care of collinear 
and infrared singularities, as well as the peaking behaviour around the W resonance. The 
sum of the soft and of the hard photon cross sections is independent of the cut-off AE. 
We have checked the independence of our numerical results from the choice of the infrared 
separator e = AE/E for 1CT 8 < e < 1CT 4 . 

The squared matrix element |-Mi| 2 has been calculated and evaluated with good agree- 
ment in several different ways, using FeynArts and FormCalc, by hand with the help of 



FORM [38], using the ALPHA algorithm |p9| , in order to check the numerical stability of the 
results. 
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3. Higher-order corrections and matching procedure 



In this section we describe the matching of the fixed EW 0(a) calculation with higher- 
order QED corrections (cfr. ref. |30]]). The latter can be included in a generic scattering 
cross section in the QED Parton Shower approach, which resums to all orders the leading 
logarithmic effects. At 0(a) the Parton Shower reproduces only the QED leading- log 
approximation of the exact EW 0(a) calculation, presented in section ^. We would like 
to combine the exact 0(a) results and QED higher orders, to improve the approximation 
intrinsic to the Parton Shower, avoiding at the same time double countings at 0(a). 

The matching of the two calculations is a non trivial task and we present, for the 
sake of clarity, the main ideas in a simplified unphysical toy model, namely in the case 
of a scattering process where only one external particle can radiate. A general expression 
of its cross section with the emission of an arbitrary number of photons, in leading- log 
approximation, can be cast in the following way: 



oo _. 

da™ = n(Q 2 , e) V - \M n , LL \ 2 d<f> n (3.1) 
z — ' n! 

n=0 

where Tl(Q 2 ,e) is the Sudakov form- factor that accounts for soft-photon (up to e) and 
virtual emissions and Q 2 is related to the energy scale of the hard scattering process. 
\-M. u ,ll\ 2 is the squared amplitude in leading- log approximation describing the process 
with the emission of n hard photons, with energy larger than e in units of the radiating 
particle energy. d§ n is the exact phase-space element of the process, with n photons in the 
final state, divided by the incoming flux factor. The cross section da°° is independent of 
the infrared separator e. 

According to the theorems for the factorization of soft and collinear photon emission, 
the leading part of the squared amplitudes can be written in a factorized form; for example 



the one-photon emission squared amplitude reads [4C] 



where 1 — z is the fraction of energy carried by the photon, I(k) is a function which describes 
the angular spectrum of the photon and P(z) = (1 + z 2 )/(l — z) is the Altarelli-Parisi 
e —> e + 7 splitting function. In eq. ( |3.2| ) we observe the factorization of the Born squared 
amplitude and that the emission factor can be iterated to all orders, giving \M n ,LL\ 2 ■ The 
Sudakov form factor H(Q 2 ,e) can be expressed as 

n(Q 2 ,e) = exp (~ 1+ log^D , 1+ = [ * dzP(z) (3.3) 



2tt m 



o 



The function I{k) has the property that J dfL, I(k) = logQ 2 /m 2 and allows the cancella- 
tion of the infrared logarithms. 

It would be desirable to include in eq. ( jO| ) the missing 0(a) contributions. The 
matching procedure can be better understood by comparing the exact 0(a) cross section 
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with the 0(a) expansion of eq. (|3.l|) , which, as we already mentioned, does not coincide, 
by definition, with an exact 0{a) result. In fact 



da. 



a,LL 



1 a r i & 
1 - — 1+ log — 

2tt m' 



|7W | 2 ^o + |-Mi,ll| 2 ^i 



= [1 + C a , LL ] \M \ 2 d$ + |.Mi,ll| 2 ^i (3.4) 

whereas an exact NLO cross section can be always cast in the form 

da a = [1 + C a ] \M \ 2 d<S> + |Mi| 2 d$i (3.5) 

The coefficient C a contains the complete virtual 0(a) and the O{o) soft-bremsstrahlung 
squared matrix elements, in units of the Born squared amplitude and |.Mi| 2 is the exact 
squared matrix element with the emission of one hard photon. We remark that C a ,LL and 
I-M^llI 2 have the same singular logarithmic structure of C a and of |.Mi| 2 . 
We observe that, by introducing the correction factors 

Fsv = l + (C a -C a>LL ), F H = i + J^lj!-J^g£ , (3.6) 



\M 



1,LL 



the exact 0(a) cross section can be expressed, up to terms of C(a 2 ), in terms of its 
leading-log approximation as 

da a = F SV (l + C a , LL )\M \ 2 d^o + F H \Mi,LL\ 2 d$i (3.7) 

Driven by eq. we can improve eq. ([O]) by writing the resummed cross section as 

CO 1 / 71 \ 

da°° = F sv n(Q 2 ,e) ^ - IJi^.i \M n , LL \ 2 d<5> n (3.8) 

n=0 ' \i=0 / 



The correction factors Fjj^ follow from the definition eq. ( |3.6| ) for each photon emis- 
sion. The expansion at 0{a) of eq. ( |3.8| ) coincides now with the exact NLO cross section 
of eq. (|3.5|). Furthermore, all higher-order leading-log contributions are the same as in 
eq. (|3.1|). It is worth noticing that Fsv,Fh,i are, by construction, infrared safe and free of 
collinear logarithms. 



Alternatively, we remark that one could have improved eq. (3.1) by adding the O(a) 
contributions missing in the leading-log approximation. However, we prefer the factorized 
formulation eq. ( p.8[ ) for different reasons. The first is that it allows to include a large 
class of radiative corrections to all orders beyond the leading logarithmic approximation. 
In fact, the factorization of the soft and collinear QED corrections is a universal property 
valid for any hard scattering process; the latter can be computed not only at tree level, 
but also including infrared safe radiative corrections; this is the case in eq. ( |3.8[ ), with the 
factors Fsv an d Fjj^. As a consequence, not only the tower of leading terms of the form 
a n log n (s/M 2 ) but also all the terms of the form a n log n_1 (s/M 2 ) log(s/A£ 2 ) (M is the 
mass of the radiating particle) are correctly taken into account by eq. ( |3.8[) . Another reason 
for choosing a factorized formulation is related to the request that the exponentiated cross 
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section should go to zero in the limit of vanishing infrared separator e: in the addittive 
formulation the missing infrared finite 0(a) terms would spoil this property. 

Equation ( |3.§| ) is our master formula for the matching between the exact 0(a) EW 
calculation and the QED resummed Parton Shower cross section. Its extension to the 
realistic case, where every charged particle radiates photons, is almost straightforward. 

It is useful to present, in the general case, the expression of the function I(k), which 
describes the leading behaviour of the angular spectrum of the emitted photons. 



i(*0 = £ QiQiwu p . P i' P J.. k E " ( 3 - 9 ) 



where Qi and pi are the electric charge fraction and the momentum of the external fermion 
I, r]i is a charge factor equal to +1 for incoming particles or outgoing antiparticles and 
equal to -1 for incoming antiparticles or outgoing particles, k is the photon momentum and 
the sum runs over all the external fermions. 
In our case the Sudakov form factor is 

U( Pi ;e) = exp(-|-/ + J dU 7 I(k)\ (3.10) 
The result of the analytical integration of I{k) is: 

f c^(A0 = -££w&ttQ;|£r5 iog4-£Q? (3-11) 

j>i i 1 • 1 j 1 1 i 

where 

Pi — PijPi 

The parameters are obtained requiring that (pi — pj) 2 = and $ —p® = PijP® —p® > 0. 
4. Hadron-level cross section 

In this section we discuss the hadron-level cross section a(pp — » lui~\-nj + X) and the pro- 
cedure to subtract the initial-state mass singularities in the calculation of the cross section 
first with 0(a) and then with the 0(a) improved with QED higher-order corrections. This 
procedure makes the resulting hadron-level cross section independent of the (unphysical) 
value of the initial-state quark masses. 

The radiative corrections are enhanced by large collinear logarithmic (mass singu- 
larities) factors, of the form a/(2n) log(s/m^). The initial-state collinear logarithms are 
universal, i.e. are independent of the hard scattering process, must be factorized out at 
an energy scale M and can be reabsorbed in the definition of the parton densities which 
describe the partonic content of the proton, like in all NLO QCD calculations. Since the 
partonic cross section described in section [2] still contains QED initial-state mass singulari- 
ties, it is mandatory to implement a subtraction procedure to avoid a double counting when 
convoluting with the proton parton densities. Recently a new set of PDFs (MRST2004QED) 
including QED effects has been published [41]. These parton densities are evolved via the 
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DGLAP equations including also the QED splitting functions. As discussed in the litera- 
ture ||, |U , the effect of the QED evolut ion on the PDFs is at the per mille level, for x 
and M 2 values probed by Drell-Yan dynamics at the Tevatron and at the LHC. 

The subtraction at 0(a) has been discussed for example in ref. [13] and is obtained 
by a redefinition of the parton densities. The hadron- level cross section at 0(a) can be 
written 

da(pp -»■ Ivi + X) = V / dxidx 2 q a {x u M 2 )q b (x 2 , M 2 ) [da + da a ] - (4.1) 

a,b J ° 

- (Aq a {x 1 ,M 2 )q b (x 2 ,M 2 )+q a {x 1 ,M 2 )Aq b (x 2 ,M 2 )) da 

where a, b run over all parton species described by the densities qi(x,M 2 ), M is the 
factorization scale, dao and da a are the Born and 0(a) partonic cross sections. The 
0(a) subtraction term is 



A qi (x,M 2 



-,M 2 ) —Q 2 

Z / Z7T 



(4.2) 

Qi and m\ are the electric charge fraction and the mass of the quark i and the function 
f(z) p3] allows to change the subtraction scheme (e.g. DIS or MS). Given the presence 



in the hadron-level cross section eq. ( |4.1| ) of the product of two parton densities, the 
subtraction procedure in a factorized form could yield terms of 0(a 2 ) which have been 
discarded for consistency at 0(a). 

It is mandatory to generalize the independence from the value of the quark masses of 
the O(a) cross section of eq. ( |4.l| ) to the cross section including also QED higher-order 



corrections. The exponentiation in the first line of eq. ([4.1D of [dao + da a ], according to 
the matching algorithm described in section |3[ would develop also higher-order initial-state 
mass singularities, which appear in the form a n log n (s/m 2 ) (n > 1). In order to remove 
the latter in a systematic way, we prefer to rewrite the partonic cross section in eq. (|4.1|); 
we split it in two terms: one free of initial-state mass singularities, da a , which will be then 
improved with the resummation, and one which contains the singular part at 0(a) , da^ b . 
We rewrite eq. (WJh by adding and subtracting the same quantity, namely 



dan + da 



da + (da a - da s a ub ) + da s a ut> = da + da, 



sub 



dat 



(4.3) 



The subtraction term is further split into a soft+virtual and a hard photon contributions: 

da s a ub = da^ v ' sub + da^ sub (4.4) 



da 



SV,sub 



da 



a 
8ira 



M 2 



mi 



Q' a Mog-y-l +Qi Mog-^-1 



M 2 



mt 



1 + z 2 



\M (s')\ 2 I sub (k) d*i = \M 



E 2 z(l-z) 1 



i=a,b 



p iC M 2 (l-(3 iC ) 2 



I+\M { z d$o= C s a ub \Mo\'d$o 



2mj 



l,sub\ 
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where the subtraction cross sections are defined in analogy to the leading log cross sec- 
tions of eq. ( ft.4|) and I su b(k) describes the initial-state radiation. In eq. ( |4.4j ) |_Mo(s')| 2 
is the Born squared amplitude evaluated at a reduced center-of-mass energy s', Pi = 
a/i — Amf/M 2 and c is the cosine of the angle of the photon with the beam axis in the 

partonic center-of-mass frame. Once integrated over the photonic variables, daa' sub devel- 
ops the same infrared logarithmic structure as da a V,sub and the sum is independent of the 
value chosen for the infrared separator. 

By applying the matching algorithm to [daQ+da a ], which is free of collinear initial-state 
logarithms, we obtain an improved hadron-level cross section including QED higher-order 
corrections. Starting from the 0(a) expression, identical to eq. (|4.1|), 



(4.5) 



dahad = / dxidx 2 q a {xi,M 2 )q b {x 2 ,M 2 )[doQ + da a ] + 

a,b J ° 

'Aq a ( Xl ,M 2 ) Aq b (x 2 ,M 2 ] 



+ q a (xi,M z )q b (x 2 ,M 2 



da 



/a (xi,M 2 ) %(x 2 ,M 2 ) 
we obtain the hadron-level QED resummed cross section 



da. 



had 



I dx x 

0,6 J ° 



dx 2 q a ( Xl ,M 2 )q b {x 2l M 2 ) x 



oo - / n \ 

%n«V) J^^y n^.M \Mu,ll\ 2 d$ n + 

71=0 n ' \i=0 / 



+ 



da sub_ { ^a^M 2 ) + Aq b{ X 2 ,M JU(jQ 



t (x!,M 2 ) %(x 2 ,M 2 ) 



dan 



(4.6) 



The variables with a tilde represent quantities subtracted of the initial-state singularities; 
more precisely the matching of section || is built using the following subtracted quantities: 



l-Mil 2 = l-Mil 2 - |-M liS „ 6 | 



\Mi,ll\ 



r~i f~iSub 



C a ,LL 



a, 



.LL 



-a 



a > 



\M^ LL \ 2 - \M^ sub \ 2 

I(k) = I(k) - I sub (k) 



The 0(a) expansion of eq. ( |4.6| ) coincides with eq. flOj ) or equivalently to eq. (|4.1| ). The 
second line in eq. ( [4.6[) is by construction free of initial-state mass singularities and we 
remark that the same property holds also for the finite 0(a) correction of the last line. 

Equation (|4.6| ) is our master formula for the computation of the hadron-level cross 
sections and the event simulation. 



5. Numerical results 

All the numerical results have been obtained using the following values for the input pa- 
rameters: 
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a = 


1/137.03599911 




= 1.16637 10 GeV 


m z = 


91.1876 GeV 


m w 


on a o r (~\ ~\ t 

= 80.425 GeV 


1 w 


= 2.124: GeV 


m H = 


i i r f A ~\ j 

115 GeV 


m e 


= 510.99892 KeV 




= 105.658369 MeV 


m T = 


1.77699 GeV 


m u 


= 66 MeV 


m c 


= 1.2 GeV 


m t = 


178 GeV 


m d 


= 66 MeV 


m s 


= 150 MeV 


m b = 


4.3 MeV 


Vud 


= 0.975 


V us 


= 0.222 


v ub = 





Vcd 


= 0.222 


Vcs 


= 0.975 


v cb = 





Vtd 


= 


Vts 


= 


v tb = 


1 



and have been computed in the input scheme described in Section & The set of parton 



density functions used to compute all the hadron-level cross sections is MRST2004QED [41] *. 
In this set of PDFs the QCD and the QED factorization scales are set to be equal and, as 
usually done in the literature |l^, [H| ], we choose M = m w , if not stated otherwise. The 
use of the PDFs set MRST2004QED implies that our numerical results correspond to the DIS 
factorization scheme. The computation of the hadron-level results requires the numerical 



evaluation of the subtraction term defined in eq. 4.2; a grid of values in the variable x, 
which is then interpolated, is obtained by means of the numerical library CUBA [42|. All the 
hadron-level results refer to the LHC, at a nominal center-of-mass energy y/s = 14 TeV, if 
not stated otherwise. 



The fine structure constant a is used instead of a^ e of eq. (^6) in the computation of 
all the 0(a) and higher-order corrections, in order to describe the real photon emission with 
the proper coupling. This replacement is formally justified in the fixed 0(a) calculation 
because it differs at 0(a 2 ). The input scheme has been implemented by computing 
all the contributions to the cross section with a = a(0), subtracting the Ar contribution, 
evaluated as well with a(0), and then rescaling the total result by (ag ee /a(0)) 2 . 

The following cuts have been imposed to select the events: 

pj_ t e > 25 GeV, p± 1 missing > 25 GeV, \r]e\ < 2.5 (5.1) 

where p±/ and rje are the transverse momentum and the pseudo-rapidity of the charged 
lepton and p± tm i ss ing is the missing transverse momentum, which in our case coincides with 
the neutrino px_ ■ 

Our results are obtained with bare muons and with "recombined" electrons. We assume 
perfect isolation of photons from the muon, which is experimentally achievable with good 
accuracy: the resulting cross sections are therefore enhanced by large muon mass collinear 
logarithms, because the photon emission is not treated inclusively in the region about 
the muon. In the case of electrons, it is experimentally difficult to separate them from the 
photon track, when the latter lyes within a cone around the lepton. We adopt the following 
recombination procedure: 

• photons with a rapidity \rj^\ > 2.5 are never recombined to the electron; 



* Since this set of PDFs includes a photon distribution function resulting from the QED evolution of the 
PDFs, the Drell-Yan cross section receives a new type of correction from 2 — > 3 photon-induced processes, 
such as, for example, 75 —> q'lvi. These contributions, which have been evaluated by Dittmaier and Kramer 
in ref. [Eq], are not considered in the present study. 
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1. 

2. 
3. 
4. 
5. 



lowest order (Born) 

final state 0(a) LL QED Parton Shower 

final state exponentiated LL QED Parton Shower 



exact 0(a) EW of eq. (4.1) 



exact 0(a) EW matched with higher-order QED corrections (best) of eq. (|4.6|) 



Table 1: Different approximations for the cross section of the Drell-Yan charged current process. 



e 


0(a) 


best 


5 • 10- 4 


4410.98 ±0.20 


4412.14 ±0.26 


i • i(r 6 


4410.84 ± 0.48 


4413.66 ± 0.9 



Table 2: Independence of the hadron-level cross section (pb) of the value chosen for the infrared 
separator s. 



• if the photon rapidity is |n 7 | < 2.5 and R ei = J yq e — r]j) 2 + < 0.1 (0 e7 is the 
angle between the photon and the electron in the transverse plane), then the photon 
is recombined with the electron, i.e. the momenta of the two particles are added and 
associated with the momentum of the electron; 

• the resulting momenta should satisfy the cuts of eq. (|5.1| ). 

In order to study the different effects of the radiative corrections on the relevant ob- 
servables, we will distinguish the approximations described in table [l]. The approximations 
2. and 3. refer to the old version of HORACE pQ ], which included the final-state QED 
corrections in leading- log accuracy (both at 0(a) and with higher orders) in a pure Parton 
Shower approach. 

5.1 Technical checks 

In this section we discuss some technical checks satisfied by the new version of HORACE, 
namely the independence from the e parameter and from the quark masses. As reference, 
we also provide the (un-subtracted) cross section at parton level for three different center- 
of-mass energies. 

The definition of the soft- and hard-bremsstrahlung corrections requires the intro- 
duction of an infrared separator e. The physical cross section does not depend on this 
parameter as shown in table ||. 

The partonic cross section is computed using quark masses to regularize the initial- 
state collinear singularities. The latter are subtracted according to the procedure described 
in Section || when calculating the hadron-level cross section, which has to be independent, 
up to terms of order m^/m\ r , of the value chosen for the quark masses. This property is 
demonstrated in table ||, both for the 0(a) and matched cross section. 
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0(a) 


best 


m q 


4410.98 ± 0.20 


4412.14 ±0.26 


m q /10 


4410.92 ±0.26 


4411.89 ±0.33 


m q /WO 


4410.99 ± 0.29 


4411.92 ±0.50 



Table 3: Indipendence of the hadron-level O(o) and exponentiated cross sections (pb) of the value 
chosen for the quark masses. 



>/i (GeV) 


60 


80 


100 


Born 


7.807 ±0.001 


5223.0 ±0.5 


25.452 ± 0.002 


0(a) 


7.274 ± 0.002 


4817.1 ±0.8 


34.476 ± 0.005 


best 


7.288 ± 0.002 


4830.6 ± 0.6 


34.208 ± 0.008 



Table 4: Unsubtracted partonic cross section (pb) at different partonic center-of-mass energies. 

Finally, the partonic cross section has been computed, within the cuts of eq. ([O]), 
without subtraction of the initial-state collinear divergences; the results are presented in 
table |. 

5.2 W transverse mass distribution 

We present in this section the numerical results for the distribution of W transverse mass, 
defined as 

M x = yfeplj P.l,v (1 - cos (f>e tU ) (5.2) 

where (f>i v is the angle between the lepton and the neutrino in the transverse plane. We 
discuss in detail the effect of the different classes of radiative corrections. Their impact, 
relative to the Born approximation, is in general of several per mille and in some cases of 
few per cent and can play a significant role for instance in the precise determination of the 
mass of the W boson m w , foreseen at the LHC with an accuracy of 15 MeV. 

In table | we show the cross section integrated within the cuts eq. (|5.l| ) and a futher 
cut on the minimum W transverse mass, varying from 50 to 2000 GeV. The relative effect 
of the 0(a) corrections grows with the cut on the minimum transverse mass, because of 
the increasing importance of the EW Sudakov logs. 

The Born results coincide for muons and electrons, up to negligible mass effects. The 
radiative corrections instead differ because of the final-state collinear logarithmic enhance- 
ment, which are absent in the case of photons recombined with the electron. All the QED 
higher-order corrections do not modify significantly the 0{a) corrections. 

In figures from ^ to ^ we show the transverse mass distribution and disentangle the 
different contributions due to the radiative corrections. In figure ^ the transverse mass 
distribution is plotted, in the range 50 < M_l < 100 GeV. The transverse mass distribution 
provides physical information in different ranges: the position of the jacobian peak and 
the shape of the distribution about the peak can be used to extract the value of the W 
boson mass, the shape of the tail of the distribution above the peak, 80 < M_l < 100 
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DUI 11 I JJ U 1 


^ + (%) 

"a \/0) 




°a \ /0 ) 


°oo \ /0 ) 


50 


t:iJ<JU.uOI I I 


-9 8 


-9 7 


-1 7 


-1 8 


100 


27.642(1) 


-5.0 


-4.9 


-3.4 


-3.4 


200 


1.79275(5) 


-7.9 


-7.7 


-6.3 


-6.3 


500 


0.084809(2) 


-14.3 


-13.8 


-12.2 


-12.2 


1000 


0.0065320(2) 


-21.9 


-21.1 


-19.4 


-19.1 


2000 


0.000273686(8) 


-32.1 


-30.5 


-28.7 


-28.1 



Table 5: Lowest-order hadron-level cross section, integrated imposing a cut on the minimum 
transverse mass and relative effects, with respect to the Born cross section, in the O(a) (8 l a ) and 
in the best (5^) approximations. 




50 55 60 65 70 75 80 85 90 95 100 

M± (GeV) 



Figure 3: Transverse mass distribution in Born, O(a) and best approximations. 

GeV, can be used to measure the W boson decay width and the large transverse mass 
tail, 200 < Mj_ < 1000 GeV, of the distribution can be an important background to the 
searches of new heavy gauge bosons. 

In figure [| we plot, in the range 50 < M± < 100 GeV the effect of the exact 
0(a) radiative correction, relative to the Born cross section, in the case of muons and 
of recombined electrons. The 0(a) contribution gives a large correction, up to ~ —10%, 
which distorts, about the W resonance, the transverse mass distribution and is responsible 
for the bulk of the shift in the extraction of the W boson mass. As shown in figure ||, in the 
range 100 < M± < 1000 GeV the EW Sudakov logarithms make the effect of the radiative 
corrections large and negative, reaching the 20% level. 

In figure ^ we disentangle, among the 0(a) contributions, the effect of all the correc- 
tions which can not be classified as QED final state-like leading-log radiation, by taking 
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^ -4 

-6 

-8 
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50 60 70 80 90 100 

Ml (GeV) 

Figure 4: Relative corrections with respect to the Born cross section due to the exact 
0(a) corrections for muons and recombined electrons final states. 




100 200 300 400 500 600 700 800 900 1000 

Mj_ (GeV) 



Figure 5: Relative corrections with respect to the Born cross section due to the exact 
0(a) corrections for muons and recombined electrons final states. 

(blue line) the difference between approximations 4. and 2. (and between 5. and 3., red 
line) of table [j] in units of the differential Born cross section. We present only the results 
for muons, being the effect similar in the electron case. We observe that they are quite 
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Figure 6: Relative effect, in Born units, of the difference between the approximations 4. and 2. of 
table 1 (blue line) and between 5. and 3. (red line). 




Mi (GeV) Mi (GeV) 



Figure 7: Relative effect on the transverse mass distribution, in Born units, of higher-order QED 
final state-like and full QED parton shower corrections. 

flat, small and negative, for M± < 80 GeV; they become larger in size and always negative 
for increasing values of M±, because of the presence of the EW Sudakov logs. From a 
comparison of figures |5] and ^, the non-factorizable weak contributions account for more 
than half of the 0(a) radiative corrections, for M± > 200 GeV. 

In figure |7| we present the effect of the higher-order (beyond 0(a)) corrections, and 
disentangle the effect of all the terms which can not be classified as QED final state-like 
leading log radiation, by considering the difference of the 3. and 2. (red line) and of 5. 
and 4. (blue line) approximations, in units of the lowest-order differential cross section. 
We present only the results for muons, being the higher-order corrections smaller in the 
electron case because of the recombination. The red line describes the effect of purely 
photonic final-state like leading log corrections, whereas the blue line represents the higher- 
order contributions of the matched cross section of eq. ( |4.6| ). The latter includes, besides 
the content of the red line, the remnant of the initial-state radiation after the subtraction 
of the initial-state singularities and the product of purely weak corrections (the Fgv factor 
of eq. (|4.6| )) with photonic radiation. Around the peak the two lines almost coincide, 
while for large M± we observe the effect of the product of the EW Sudakov logs times the 
0(a) photonic correction. 
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scheme 


Born 


0(a) 


f (%) 


a(0) 


4244.68 ± 0.09 


4360.5 ±0.6 


±2.73 




4536.03 ± 0.07 


4411.0 ±0.2 


-2.76 



Table 6: Born and 0(a) hadron-level cross sections (pb)and effect of the 0(a) corrections, ex- 
pressed in units of the corresponding Born cross section, in the a(0) and in the G M schemes. 

6 
4 
2 


g " 2 

<« -4 
-6 
-8 
-10 
-12 

50 55 60 65 70 75 80 85 90 95 

M± (GeV) 

Figure 8: Relative corrections to the transverse mass distribution in the and in the a(0) 
schemes, expressed in units of the corresponding Born distributions, and their difference. 

As we already discussed in Section H we can compute the cross sections in the or 
the a(0) input scheme. In table ^, we compare the cross sections obtained in the two input 
schemes, in Born and in 0(a) approximations and the corresponding relative corrections. 
The difference between the cross sections in the two schemes is reduced when going from 
the Born to the 0(a) approximation and amounts to about 6% (Born) and 1% (0(a)), 
respectively. The relative correction in the two schemes is of the same order (~ 3%) but of 
opposite sign. This can be understood taking into account that, as previously discussed, 
in the scheme, at a variance with the a(0) scheme, universal virtual corrections are 
absorbed in the lowest-order cross section. It is worth noticing that the 0(a) corrected 
transverse mass distribution differs in the two input schemes as shown in figure |8|, where 
we plot the relative corrections in the two schemes in units of the corresponding Born 
distributions and their difference. 

Another source of uncertainty, which is not of purely EW origin, is the choice in the 
parton densities of the factorization scale M. In order to study this dependence, we set 
M = £,m w and consider the canonical range 1/2 < £ < 2. We define the two following 
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Figure 9: Relative effect, with respect to the corresponding Born, of the 0(a) corrections, com- 
puted with the factorization scale M = m w /2 and M = 2m w . 
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Figure 10: Relative effect of the 0(a) corrections, computed with the factorization scale M 
m w /2 and M — 2m w , expressed in units of the Born cross section with M = m w . 



relative corrections: 



6(M) 



<Ta(M) 

MM) 



A(M) 



<Jg(M) -<7 (M) 



(5.3) 



In figure ||| we plot, for the transverse mass distribution, 6(0. bm w ) and 5(2m w ). The 
difference between the two curves can be interpreted as mainly due to the dependence of 
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the 0(a) cross section on the choice of the QED factorization scale. We observe a variation 
at the per mille level of the transverse mass distribution, as already remarked in ref. pi 



In figure |T0] we plot, for the transverse mass distribution, A(0.5m w ) and A(2m w ). 
The difference between the two curves can be interpreted as an estimate of the uncertainty 
due to missing 0(aa s ) corrections, which are of the order of 1% around the W resonance. 
In fact the numerator in A(M) is proportional to the hadron-level 0(a) corrections, which 
are a QCD-LO result. A change of the scale in the PDFs gives an estimate of the QCD- 
NLO 0(aa s ) effects t. We plot in fi gure 10 the 0(a) correction at two different scales, 
normalized to the Born cross section evaluated at a fixed scale M = m w . In fact it is 
known that the introduction of QCD corrections stabilizes the total cross section against 
scale variations and we mimic this effect by setting M = m w in the denominator of A(M). 
In this way the uncertainty due to missing 0(aa s ) corrections can be disentangled from 
other purely QCD effects. 

To conclude this Section, we would like to comment on the possible relevance of the 
different effects and uncertainties on the W boson mass measurement. It is known in the 
literature Q [l4|, that the distortion of the M± distribution around the peak due to 
the EW 0(a) and QED higher-order final-state corrections induces a shift of the extracted 
m w of the order of 100 and 10 MeV respectively (for the muon case). In the new version 
of the generator HORACE both effects are included in a unique tool. 

We remark that the theoretical uncertainty due to the QCD factorization scale choice, 



shown in figure |10|, modifies the impact of the EW radiative corrections by an amount 
which can be of the order of 1%: these effects may induce a systematic error in the m w 
measurement which can be comparable with the aimed experimental accuracy and should 
be carefully considered in future experimental analyses. 

5.3 Rapidity distributions and charge asymmetry 



In figure 11 the muon and electron pseudo-rapidity distribution is presented in the approx- 
imations 1., 4. and 5. of table [l|. We observe that the effect of the radiative corrections is 
almost constant over the whole range in pseudo-rapidity and that it is dominated by the 
0(a) term, which gives a correction negative of approximately -2.7% in the case of muons 
and of -1.8% for recombined electrons. Higher-order terms modify the result at the per 
mille level. 

In figure |l^ the VF-boson rapidity is also presented. With the chosen cuts, this distribu- 
tion is essentially flat in the central rapidity interval \yw\ < 1-7. The radiative corrections 
are negative and quite flat and reduce the Born distribution of about -2% for the electrons 
and of -3% for the muons, as shown in the inset. 

The W charge asymmetry presented in figure [l3| is defined as 

da+/dr] e - da-/dr] e 
A{m) = da+fdnt + da-/d Vi (5 ' 4) 



+We are not considering in this analysis other contributions of the same perturbative order, like the one 
due to the QED PDFs evolution multiplied with the diagrammatic QCD NLO corrections, because they 
are numerically suppressed. 
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Figure 12: W boson rapidity distribution in Born, O(a) and best approximations. 

where da^ = da(pp — > l^v + X); the asymmetry is due to the partonic content of the 
incoming hadrons, which leads to different lepton pseudo-rapidity distributions in the pro- 
duction of W + or W~ . The charge asymmetry can be studied both at the Tevatron and at 
the LHC, with different results due to the two different initial states and to the different 
ranges of partonic x probed at the two colliders. At the LHC the function is even under 
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Figure 13: Charge asymmetry as function of the muon pseudo-rapidity at Tevatron (y/s = 1.96 
TeV) (left panel) and LHC (y/s = 14 TeV) (right panel), in Born and O(a) approximation. 
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Figure 14: Distribution of the W boson transverse momentum, defined as the transverse momen- 
tum of the final state lepton pair. 

rji — ► — rjt, whereas it is odd at the Tevatron. The effect of the 0(a) corrections is at the 
1% level, while higher-order effects are numerically negligible. 

5.4 W transverse momentum and photonic observables 

Real photon radiation gives to the final state lepton pair a transverse momentum, which 



defines the W boson transverse momentum, whose distribution is presented in figure 14 



in the approximations 2., 3., 4. and 5. of table [T] The large difference in the tail is 



*In the present study, the transverse motion of the W boson, as due to initial-state QCD radiation, 
is neglected, because it requires a careful inclusion of QCD corrections, which is beyond the scope of the 
paper. 
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Figure 15: Distribution of the hardest photon transverse momentum. 



due to a better description of the hard photon radiation given by the exact 0(a) matrix 
element, with respect to its Parton Shower approximation. The same comment applies 
when including, in the two cases, multiple-photon radiation. 

We now present some distributions for the radiative Drell-Yan event, with a hard 
photon associated to the large transverse momentum lepton pair. This signature can 
be a source of information to study for instance the trilinear gauge boson WWj vertex 
1 45]. We select the events by imposing the cuts of eq. (5.1) and requiring that the most 
energetic photon is detected, i.e. with rapidity |r/ 7 | < 2.5 and Ej > 3 GeV. The transverse 



momentum and the rapidity distributions of the hardest photon are plotted in figures 15 



and 16. As in the W transverse momentum case, the large difference in the tail in figure 15 



is due to the better description of the hard photon radiation given by the exact O{o) matrix 
element. Concerning the hardest photon rapidity distribution in figure 16, we compare first 
the lowest order results, which have been obtained in two different approximations, namely 
2. and 3. of table |l| and observe the increase of the cross section given by the exact matrix 
element. When including multiple-photon radiation, the lowest order results receive a 
negative correction. 



6. Conclusions 



In this paper we presented a precision calculation of the charged current Drell-Yan process, 
which includes both the exact EW 0{a) and the leading-log multiple-photon corrections. 
In order to keep consistently under control the two effects, a matching algorithm between 
the fixed order calculation and a QED Parton Shower, avoiding the double counting of 
the leading-log corrections, has been devised and implemented. To our knowledge, the 
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Figure 16: Distribution of the hardest photon pseudorapidity. 



matching algorithm here presented is the first example of such an application in the field 
of EW radiative corrections. Initial-state QED collinear singularities have been regularized 
by means of finite quark masses, requiring a subtraction of the initial-state logarithms (in 
analogy to QCD NLO calculations) which are already accounted for in the evolution of 
the PDFs. The subtraction procedure, already known at 0(a), has been generalized to 
the QED resummed cross section, to systematically remove initial-state logarithms at all 
orders. Our results have been verified to be completely independent of the value of the 
quark masses, as expected. 

We studied, with the new version of the event generator HORACE, the impact of different 
classes of radiative corrections on several physical observables, showing the importance of 
combining in a unique tool fixed order results with the resummation of multiple-photon 
radiation. In fact, radiative corrections induce effects ranging from several per mille to few 
per cent of the lowest-order cross section; they change the integrated cross section and the 
shape of the distributions. 

For example, the value of the W boson mass, extracted from the transverse mass 
distribution, is shifted, in the case of a final-state muon, by EW 0(a) corrections of 
0(100 MeV) and of O(10 MeV) by multiple photon emission. These effects are important in 
view of the measurement foreseen at the LHC and Tevatron with an accuracy of Am w ~ 
15 and ~ 30 MeV and can now be included in a systematic way in the analysis of the 
experimental data. Far from the Jacobian peak, the tail of the transverse mass distribution, 
where the Drell-Yan process is a background to new gauge boson searches, receives large 
negative corrections due to the presence of the 0{a) EW Sudakov logarithms, of the order 
of 15-30%. The effect of EW corrections is important also for the estimate of the detector 
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acceptances in view of exploiting the W production process as a precise luminosity monitor 
at the LHC, with an accuracy of some per cent. 

We also studied radiative events, where at least one hard photon accompanies the 
final state lepton pair, which are a useful source of information for instance to measure 
the trilinear gauge boson vertex. In view of precision studies, the effects of the radiative 
corrections and of treating with the exact matrix element the photon radiation are sizeable 
on observables like the photon transverse momentum and rapidity distributions. 

Besides the EW 0(a) and multiple-photon corrections, we studied some of the re- 
maining theoretical uncertainties. We considered the ones due to the EW input scheme 
choice and to the QCD factorization scale choice. We would like to remark that they 
should be considered in detail to quantify, for example, the systematic error on the W 
mass measurement. 

A possible development of the work presented in this paper is the implementation in 
HORACE of the complete EW 0(a) corrections to the Z production process and a detailed 
phenomenological study of their impact on the neutral current Drell-Yan process, which is 
important for many measurements and calibrations at hadron colliders. Furthermore, we 
are now working to combine the results presented here with QCD corrections, aiming at 
providing a unified tool which includes the relevant EW and QCD effects and which can 
be very useful for the experimental collaborations. 
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